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ABSTRACT 

We have retrieved Spitzer archive data of pre-stellar cores taken with the Multiband 
Imaging Photometer for Spitzer (MIPS) at a wavelength of 160 /im. Seventeen images, 
containing eighteen cores, were constructed. Flux densities were measured for each 
core, and background estimates were made. Mean off-source backgrounds were found 
to be 48±10 MJy/sr in Taurus and 140±55 MJy/sr in Ophiuchus. Consistency was 
found between the MIPS 170/xm and ISOPHOT 160/im calibration. Fourteen cores 
were detected both by MIPS and our previous submillimetre surveys. Spectral Energy 
Distributions (SEDs) were made for each core, using additional 24- and 70-^m data 
from the Spitzer data archive, as well as previous infra-red and submillimetre data. 
Previous temperature estimates were refined, and new temperature estimates were 
made where no Infrared Space Observatory {ISO) data exist. A temperature range of 
8-18K was found for the cores, with most lying in the range 10-13K. We discount 
recent claims that a large number of pre-stellar cores may have been misclassified and 
in fact contain low luminosity protostars detectable only by Spitzer. We find no new 
protostars in our sample other than that previously reported in L1521F. It is shown 
that this has a negligible effect on pre-stellar lifetime estimates. 
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1 INTRODUCTION 

One of the greatest challenges currently facing star forma- 
tion is the origin and nature of the cold dense molecu- 
lar cloud cores where protostars form. Interstellar dust in 
these cores absorbs radiation at ultraviolet and visible wave- 
lengths and re-radiates the energy in the far-infrared and 
submillimetre. This series of papers has used the great sub- 
millimetre and infrared telescopes to trace that emission and 
to explore the initial conditions of the star formation pro- 
cess.^ 

iBeichman et al.l l|l986l ) compared the catalogue of dense 
cores surve yed in molecular lines b y Myers and his collabo- 
rators (see iBenson fc Mveri 1 19891 . and references therein) 
with the IRAS point source catalogue and identified a 
new class of dense molecular core. The lack of an in- 
frared point source is taken as evidence that there is 
no embe dded pr o tostellar heat source in the se 'starless' 
cores. In iPaper j l|Ward- Thompson et al.l Il994h of this se- 
ries we termed those starless cores which we believed 
to be undergoing star formation 'pre-protostellar' cores 
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(subsequ entl y shortened to 'pre-stellar' for brevity). In 
(Paperl? and 'Paper llf ('Andre. Ward- Thomps on, fc Mottd 
1996; Ward-Thompson, Motto, & Andre 1999) we probed 
the 1.3m continuum emission from nine pre-stellar cores 
with the TRAM telescope. 



In lPaper iv| ||jessop fc Ward-Thompsonll200ll ) a radia- 
tive transfer model fitted to C^^O observations of the core 
L1689B showed that there was freeze-out of CO towards the 
centre of the core, and that the combination of molecular 
line and continuum observations could constrain the tem- 
perature and density profiles of pre-stellar cores. The freeze 
out of molecules with density is widely established with, for 
example, similar results being found in the dense cores L1544 
and B68 l|Caselh et al.ll2002inBe"rgin et al.ll2006l ). 



The submillimetre and millimetre spectral energy dis- 
tribution (SED) of the observed pre-stellar cores from our 
first four papers showed that our observations were on the 
Rayleigh- Jeans side of a greybody spectrum that peaks 
somewhere around 100-200/im. Low atmospheric transmis- 
sion means that observations of this peak are normally im- 
possible fro m the ground (with the nota ble exception of 
THUMPER, IWard-Thompson et all |2005| ). and observers 
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are forced t o fall back on a irborne (e.g. Kuiper Airborne 
Observatory. ICameronllToT^ or space based observatories. 

The Infrared Space Observatory {ISO) was a 65- 
cm diameter space te lescope designed t o observe be- 
tween 2.5 and 240/im (jKessler et al. Il996l ). It was oper- 
ated by the European Space Agency (ESA) and ceased 
operations in early 1998. We used the imaging photo- 
polarimeter (ISOPHOT) on ISO to search for the peak 
of the pre-stellar greybody spectrum by imaging eighteen 
cores at 90, 170 and 200/jm and presented our results in 
iPaper Vl (|Ward- Thompson. Andre, fc KirkI 120021 '). For each 
core maps of colour temperature were constructed using the 
pixel-to-pixel 170/200-/im spectral index. These showed that 
most of the cores are either isothermal or are warmer at the 
edges, consistent with them lacking an internal heat source 
(i.e. a protostar) and with their surfaces being heated by the 
local interstellar radiation field. Full radiative t ransfer anal- 
ysis also shows similar temperature g radients (|Evans et al.l 
l200ll : IStamatellos fc Whitworthll200a ). 

Spectral energy distributions (SEDs) were constructed 
using the ISOPHOT data and our previous millimetre and 
submillimetre data. These were used to estimate the energy 
budget of the cores and showed that their infrared luminos- 
ity was comparable to the incident energy at shorter wave- 
lengths. This confirms that the cores are indeed pre-stellar 
and are only heated by the local interstellar radiation field. 
Our 200 and 170^m ISOPHOT observations showed the im- 
portance of far-infrared data to the accurate determination 
of the s pectral en e rgy distribution of pre-stellar cores . 

In iPaper Vj l|Kirk. Ward- Thompson, fc Andrei 120051 ) 
we used the Submillimetre Common User Bolometer Array 
(SCUBA) on the James Clerk Maxwell Telescope (JCMT) to 
conduct a survey towards 52 pre-stellar cores a t 450 and 850- 
fira (c.f. lShirlev et al.ll200Gl : lYoung et al.ll2006l ). Twenty-nine 
of the cores were detected and were separated into 'bright' 
and 'intermediate' groups based on their 850-/xm peak inten- 
sities. This s plit matches the list o f evolved prestellar cores 
identified by ICrapsi et aP (|2005al) using N2H+ and N2D+ 
obervations. The ISO data from Paper vl was used to calcu- 
late the physical parameters of the cores. This showed that 
while the normalised column density profiles of the cores 
have the same shape as an unstable Bonnor-Ebert Sphere 
l|Ebertlll955l :lB onnoj [T956I ) the implied physical character- 
istics of the spheres did not match those calculated for the 
cores. 

The Multiban d Imaging Photometer for Spitzer (MIPS, 
iRieke et al. I I2OO4I ) is a far-infrared c amera on the Spitzer 
Space Telescope jWerner et al.l |2004| ) . Spitzer is an 85-cm 
diameter cryogenically-cooled satellite telescope designed to 
operate from 3 to 160 ^m. MIPS provides Spitzer's long 
wavelength capabilities and has imaging wavebands at 24, 
70, and 160 with telescope-limited resolutions of 6, 18, 
and 40 arcsec respectively. The sensitivity, resolution and 
large field of view allows for the unparallelled rapid mapping 
of extended star formation regions. 

In this paper we present MIPS 160-/im images and 24- 
& 70-^tm flux densities of 18 regions around pre-stellar cores 
that we have downloaded from the Spitzer data archive. In 
Section [2] we outline our search strategy and post-pipeline 
processing. In Section [3] we present and discuss the im- 
ages and extract flux de nsity mea surements. In Section |4] 
we reprise our SEDs from IPaper Vl and update them by in- 



cluding the new Spitzer data from this paper. In Section [5] 
we discuss the implications for the nature and lifetimes of 
pre-stellar cores. 



2 OBSERVATIONS 

2.1 Post-pipeline processing 

A search of the Spitzer d ata archive was m ade using the 
LEOPARD software tool |[LEOPARD||2005| '). The positions 
of pre-stellar cores from Paper Vll were cross-referenced 
against the archive and coincident 160-/im observations were 
identified. The unflltered basic calibrated data (BCD) ver- 
sion of the MIPS 160-/^m data were downloaded and re- 
gridded onto images with 15 arcsec diameter pixels using 
the MOPEX routine MOSAIC (Makovoz fc Khan 2005). 
Most of the archive data had been reduced with software 
pipeline version Sll while the region L1517 was reduced 
with pipeline version S13.2. 

The MIPS 160-/im array only has half the fleld-of-view 
of the larger 24 and 70- /^m arrays. For the slow and medium 
raster scanning modes this is not a signiflcant problem as 
the field-of-view is stepped forward in increments smaller 
than the width of the 160-/im array. However with the fast 
raster scanning mode, used to efficiently observe large ar- 
eas, the array is stepped forward in increments equivalent 
to the width of the shorter wavelength arrays, resulting in 
the sky plane being under-sampled by a factor of 2 at 160- 
^m. Thus the recommended best practise for fast-scanning 
is for 4 coverages over each point — an offset forward and 
backscan over the same area to fully sample the sky and a 
complete repeat observation to r eject spikes caused by cos- 
mic ray strikes on the detectors (jCordon et al.|[2005l ). 

Even with this redundancy it is necessary to reject 
outlying measurements during the re-gridding and to clean 
the fast-scanned maps afterwards. The KAPPA task FF- 
CLEAN was used iteratively to reject pixels that lay more 
than 3(T way from the local pixel mean and then the task 
FILLBAD was used to fill the holes left b y rejected pixels 
with a mean of the four- nearest good-pixels (jCurrie fc Berrvl 
l2004h . The percentage of pixels filled by this method was 
found to be a product of scan-speed and absolute back- 
ground intensity. The intensity of the good pixels was not 
altered by this method. 

The medium scan maps did not require flltering while 
the fast scan areas required less that 30 per cent filling with 
all four-pixels for the mean coming from the eight neigh- 
bouring pixels. This gives a filtering scale of 3 pixels, 45 
arcsec, comparable to the 40 arcsec MIPS beam at 160-/im. 
Filtering in this manner will change the resolution, but only 
for those pixels actually infilled. We can estimate the mag- 
nitude of this change by adding the filtering length to the 
intrinsic resolution in the normal quadrature manner, but 
we weight the filtering scale by the percentage of pixels re- 
placed. Using an upper-limit of 30 percent gives an estimate 
of the "effective resolution" of 47 arcsec, which is 3.1 pixels. 

Comparison of pixel statistics before and after filtering 
show that within the 50 per cent contour the the filtering 
did not significantly alter the mean pixel value, this is to 
be expected as the flUing was based on the local mean. The 
peak intensity within that contour was also unaffected as, by 
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Table 1. Spitzer archive details. The Data Set Identifiers are 
given for all Spitzer archive data included in this paper. 



Region 




Data Set Identifier 


L1498 




12026624 


L1521C & 


I SMM 


11229696/9952 


L1521B 




11229184/9440 


L1521A & 


z D 


11228672/8928 


L1521E & 


s F 


11228160/8416 


L1517 




12662784 


L1512 




12021248 


L1696A 




5759488, 5767680, 4321280/1536 


L1709A 




5758720, 5766912, 4321536/2048 


L1689A 




5758464/7952, 5766656/144, 4322048 


L1709C 




5757952, 5766144, 4321792 


L1689B 
B68 




5749504, 5753867, 4321792 
12023552 


B133 




12023808 



definition, tlie infill value cannot be greater tlian tiiat of any 
of the surrounding pixels. However, this does not preclude 
the true intensity peak being coincident with a blank pixel. 
The upper-limit of the number of pixels infilled is consistent 
with a 70 percent certainty of measuring the true peak in- 
tensity. The only core at variance with the above is L1689A 
which was close to the absolute saturation limit (see follow- 
ing section). Approximately 60 per cent of its area had to 
be reconstructed and it is possible that the peak intensity 
has been underestimated. This combination of FFCLEAN 
and FILLBAD was found to be relatively simple, robust 
and fast. 

Table [1] lists the 'Data Set Identifier', a unique track- 
ing number used to identify individual Spitzer observations, 
for each of our fields. The regions analysed in this pa- 
per were all taken with the MIPS fast-scanning mode, ex- 
cept for L1498, L1512, L1517, B68, and B133, which were 
taken with the medium-scanning mode. The Ophiuchus re- 
gions (L1689A/B, L1709A/C, and L1696A) are part of the 
area mapped by t he 'Cores to Disks' (C2D) legacy survey 
vans et al.ll2003l ). although the BCDs for our maps were 
downloaded and processed by us independently of the C2D 
data release of this region. However, cross-checks have shown 
that the two reconstructions are consistent. 



2.2 Detector Saturation 

The Spitzer 160- /^m detector is a 2x20 pixel stressed Ge;Ga 
photoconductor array. During an integration the rate of volt- 
age increase on a pixel is measured by non-destructively 
reading the array 8 times per second. The 'slope' of volt- 
age with time is proportional to the incident intensity and 
is fitted by the automate d software pipeline using pixel-by- 
pixel least-squares fitting (|Gordon et al.|[2005l ) . Responsivity 
changes caused by bright sources and cosmic ray strikes are 
tracked by simulator flashes every two minutes. The abso- 
lute calibration is checked against celestial standards at the 
start of each observing campaign (MIPS 2006). 

The voltage ramp from a sufficiently bright source will 
saturate the instrument electronics before the end of the in- 
tegratio n cycle . The saturation speed at 160/im is 20 MJy/sr 
in 10 s iMIPSllioO^ 'l. The fast-scanning method has a per 
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Right Ascension (2000) 



Figure 1. Map of data coverage for a part of Ophiuchus. The 
greyscale shows the relative number of measurements that were 
propagated to the final map, black is low, white is high. The 
contours show intensity from the output map. The black contour 
is at 200 MJy/sr and the white contour is at 300 MJy/sr. 

point integration time of 3 s so will begin to saturate above 
~70 MJy/sr. The number of useful readouts, and therefore 
the goodness of the slope calibration, decreases with increas- 
ing source brightness until the source instantly saturates the 
detector and makes an estimation of its brightness impossi- 
ble. This does not mean that usable data cannot be taken 
during saturation, but there must be a sufficient number of 
voltage readouts to characterise the slope. 

We can estimate the magnitude of the saturation prob- 
lem by inspecting the number of slope measurements for 
which a valid intensity calibration could be made on a point- 
by-point basis. Figure [T] shows a map of the relative num- 
ber of valid slope measurements superposed on contours of 
sky brightness. There is a distinct anti-correlation between 
source brightness and the relative number of valid measure- 
ments with the effect first becoming noticeable above 200 
MJy/sr. The effect can be compensated for by increasing 
the number of observation repeats, but absolute saturation 
will still occur. From inspection of the data this absolute 
saturation appears to occur above 1000 MJy/sr. These in- 
tensities correspond to respective saturation times of 1 and 
0.2 seconds — the equivalent of ~8 and ~2-3 readouts - 
with the brighter intensity saturating in the shorter time. 

Of the cores included in this paper (see the following 
Section) only those from the Ophiuchus region regularly ex- 
ceed 200 MJy/sr. L1709 A and C only just exceed this, while 
L1689B and L1696A double it. Only L1689A approaches the 
apparent level of absolute saturation. However, as we shall 
show in Section |4l all their flux densities are consistent with 
those measured by ISOPHOT at 170/im. Therefore we be- 
lieve the data reduction and calibration holds good for all 
of the data presented in this paper. 



3 RESULTS 

3.1 Source Positions and Flux Densities 

Table [2] lis ts the Spi tzer positions of the SCUBA core re- 
gions from iPaper Vll that we found had been observed by 
MIPS and stored in the Spitzer data archive. The seventeen 
regions that were downloaded cont ained nine teen separate 
SCUBA core search positions from iPaper Vll . The SCUBA 
non-detection of L1517D is confirmed by MIPS not to be a 
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Table 2. Spitzer source positions, background intensities, and source fiux densities. The quoted positions are the centroids of the 150- 
arcsec apertures used to measure the flux densities in Columns 7 and 8. Sources are listed in order of increasing Right Ascension. Column 
4 lists the absolute peak intensity measured near the centroid before any background correction has been applied. Columns 5 and 6 list 
the mean and standard deviation of the source backgrounds as measured near to the minimum of the map (off-cloud) and close to the 
embedded core (on-cloud). Columns 7 and 8 list the background-subtracted flux densities in a 150 arcsec aperture centred on the centroid 
position using the backgrounds shown in Columns 5 and 6. The quantities in Columns 4—8 are quot ed to 2 signi ficant figures, the errors 
are quoted to the same number of decimal places, and the absolute calibration error is 20 per cent ||MIPS||2006|) . The error on the peak 
intensity (Column 4) is equal to the larger of the errors on the off-cloud or on-cloud background and is quoted to the same significance 
as the peak intensity. The error on the two backgrounds is the pixel-to-pixel standard deviation of a fiat r egion, of s ize comparable to 
the flux density aperture, at the level of the background. Column 9 lists the SCUBA detection group from lPaper Vj — B is 'bright', I 
is 'intermediate' and n.d is 'non-detection.' 



Core 


Right Ascension 


Declination 


jP&ak 
^160 


Backj 


;round 




•S'leo 


SCUBA 






(2000) 


(2000) 


off-cloud 


on-cloud 


off-cloud 


on-cloud 


Group 












[MJy/sr] 


[MJy/sr] 


[MJy/sr] 


[Jy] 


[Jy] 




L1498 


04'' 


lO"' 


53.8^ 


+25°09'12" 


82±2 


35±2 


57±2 


16±1 


6.8±0.1 


I 


L1521C 


04'' 


19™ 


22.6* 


_|_27°14'49" 


120±10 


49±3 


83±3 


23±1 


9.2±0.1 


n.d 


L1521SMM 


04'' 


21™ 


00.4= 


-|-27°02'32" 


<81 


47±3 


69±4 


3.2±0.1 


0.86±0.10 


I 


L1521B 


04'' 


24™ 


15.9= 


-f26°37'21" 


87±3 


45±3 


64±3 


14±1 


5.9±0.2 


n.d. 


L1521A 


04'' 


26™ 


43.1= 


-|-26°16'00" 


<80 


47±3 


71±3 


11±1 


0.54±0.15 


n.d. 


L1521D 


04'' 


27™ 


48.6= 


+26°18'04" 


82±4 


46±2 


63±4 


lOitl 


3.4±0.2 


B 


L1521F 


04'' 


28™ 


37.9= 


+26°51'08" 


140±10 


59±4 


110±6 


25±1 


5.6±0.3 


B 


L1521E 


04'' 


29™ 


15.0= 


-1-26'" 14' 10" 


94±4 


39±2 


78±4 


17±1 


1.0±0.2 


I 


L1517A 


04'' 


55™ 


17.0= 


-f30°32'53" 


llOitlO 


69±3 


85±3 


15±1 


8.9±0.1 


I 


L1517B 


04'' 


55™ 


15.5= 


+30°37'38" 


lOOitlO 


69±3 


85±3 


9.8±0.1 


3.2±0.1 


B 


L1512 


05'' 


04™ 


10.7= 


-|-32°42'56" 


lOOitlO 


49±3 


68±2 


18±1 


lOitl 


I 


L1696A 


16'' 


28™ 


29.3= 


-24°16'51" 


470±40 


130±40 


290±20 


llOilO 


45±1 


B 


L1709A 


16'' 


30™ 


47.8= 


-23°41'36" 


220±10 


140±10 


ISOitlO 


26±1 


9.8±0.5 


I 


L1689A 


16'' 


32™ 


15.8= 


-25°03'40" 


910±30 


280±20 


450±30 


200±10 


130±10 


I 


L1709C 


16'' 


33™ 


59.3= 


-23°40'18" 


230±10 


95±9 


140±10 


46±1 


25±1 


n.d. 


L1689B 


16'' 


34™ 


50.6= 


-24°38'15" 


460±30 


160±20 


290±30 


88±1 


37±1 


B 


B68 


17'' 


22™ 


38.7= 


-23°50'23" 


140±10 


73±4 


96±5 


21±1 


12±1 


I 


B133 


19'' 


06™ 


09.8= 


-06°53'04" 


160±10 


63±2 


93±6 


32±1 


19±1 


B 



real core (see below), so hereafter we limit our discussion to 
the remaining 18 cores. For those 18 cores we created maps 
that matched the approximate area (^18 ar cmin square) 
of our ISOPHOT observations from IPaper vl These maps 
(see below) contain additional sources to those mapped with 
SCUBA, however, in this paper we restrict ourselves to an 
examination of the SCUBA selected sources and leave anal- 
ysis of the other objects to a later paper. 

A 150-arcsec diameter aperture was placed at each 850- 
/xm position and a centroid was performed to determine the 
equivalent 160-/^m position — these are the positions listed 
in Columns 2 and 3 of Table [51 In most cases the coordinates 
are coincident with the SCUBA positions to within a differ- 
ence equivalent to the effective resolution of the MIPS maps 
(47 arcsec, see above). The maximum pixel value in the re- 
gion of the centroid was taken to be the peak intensity of the 
core and is listed in Column 4. The background-subtracted 
(see below) flux densities measured from the centroid aper- 
tures are listed in Columns 7 and 8. The core L1521SMM has 
two peaks at 850/xm, but only one of these was detected at 
160/im. Given their close proximity we use half-sized aper- 
tures to avoid overlap and list positions determined from the 
SCUBA maps. 

A major difference between the operation of MIPS on 
Spitzer and SCUBA on the JCMT is that the JCMT has 
a chopping secondary mirror which is used to reject bright 
atmospheric emission. This creates a problem when trying 
to compare data from both telescopes as Spitzer is sensi- 
tive to the entire dynamic range and spatial scale of the 
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Figure 2. A comparison of the off-cloud and on-cloud MIPS 160- 
/jm backgrou nds with the equivalent ISOPHOT 170-/im back- 
grounds from lPaper Vl . The two fan-lines denote the ±25 per cent 
envelope that would be expected from the 30 percent ISOPHOT 
and 20 percent MIPS calibration errors. 



source emission, but SCUBA is only sensitive to the relative 
brightness and spatial scales within a region determined by 
the chop throw of the secondary mirror (see section 5.4 of 
iPaper Vll) . It is therefore necessary to subtract a background 
from each Spitzer image in order to directly compare obser- 
vations from each instrument. 

In the manner of our earlier ISOPHOT analy sis we es- 
timate background intensities from the images. In lPaper vl 
we estimated three different background levels from the 
ISOPHOT maps - background 1 (Bl) was the lowest pixel 
in the map, background 2 was estimated from the extended 
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emission that surrounded the closed contours of the core, 
and a third "best-fit" background was estimated by fitting 
a slope to the same extended emission. MIPS and ISOPHOT 
have comparable beam sizes, but ISOPHOT undersampled 
its beam as a trade ofT for increased sensitivity. The in- 
creased information in the MIPS maps makes it harder to 
fit a simple 2D plane to the background so we forego this 
method in fav our of th e other two background subtraction 
methods from iPaper Vl 

In order to duplicate the first ISOPHOT background 
we constained our map sizes to the same approximate map 
size as the ISOPHOT observations. We then measured an 
"off-cloud" background across the lowest region of the map. 
The off-cloud background and standard deviation for each 
core is listed in column 5 of Tableland is shown as a dotted 
contour on the maps in Figures [3] and |4l The left panel of 
Figure [2] shows a graph of the ISOPHOT Bl background 
against the new off-cloud background. All points except one 
outlier (L1689A) are within the ±25 per cent calibration 
error fan (based on a mean of the 30 per cent ISOPHOT 
and 20 per cent MIPS calibration errors). 

The agreement between backgrounds shows that both 
instruments have consistent calibrations, but to be able to 
compare our measured flux densities with chopped SCUBA 
data we need to subtract a greater background level. To 
allow this an "on-cloud" background was measured at the 
level of the extended or filamentary material which appears 
to surround the closed contours of each core. A further con- 
straint on this background is that it has to be comparable 
in distribution to the size of the SCUBA chop throw (2 ar- 
cmin in this case) and in this aspect it is analagous to a 4 
arcmin diameter sky annulus. The on-cloud background and 
standard deviation are listed in column 6 of Table [2] and are 
shown as the lowest solid contour on the maps in Figures 
and |4l The on-cloud background is comparable to the B2 
ISOPHOT background; the right panel of Figure [2] shows 
a graph of the on-cloud background against the B2 back- 
ground. As with the off-cloud case these two backgrounds 
are in agreement. 

The peak intensities listed in column 4 of Table [2] are 
dominated by the intensity of the background with the on- 
cloud background contributing, on average, two-thirds to the 
full range of the peak intensity. The standard deviations of 
the two backgrounds are estimated from the pixel-to-pixel 
variation across the relatively flat region used to estimate 
each background. 

The majority of the data in our sample comes from 
two distinct star formation regions — the Taurus region at 
~04 hours R.A. and the Ophiuchus region at ~16 hours 
R.A. The mean off-cloud background showed a significant 
difference between these two regions. The weighted mean off- 
cloud background was 48±10 MJy/sr i n Taurus and 140±55 
MJy/sr in Ophiuchus. From [Paper Vl the equivalent mean 
Bl backgrounds at 170-/^m were 58±8 MJy/sr in Taurus 
and 142±41 MJy/sr in Ophiuchus. This further shows that 
we are probing the same emission with each instrument and 
that their absolute calibrations are consistent. 

3.2 Imaging Data 

Figures |3] and |4] show the 18 cores that were downloaded 
from the Spitzer data archive. The greyscale shows the 160- 



Hm dust continuum and is scaled relative to the dynamic 
range of each field to avoid domination from strong sources. 
The dark contours plot the same continuum. The lowest 
solid contour is at the level as the on-cloud background (see 
Table [2] and above) , the other contours are then separated 
from it by increments of 2-a (solid above it, dashed below). 
The lowest contour is a dotted contour that shows the level 
of the off-cloud background. Superposed on each map are 
white contours that show the sm oothed 850-^m S CUBA 
continuum data from lPaper v3 and lKirk et al] l|2006l ) , which 
cover a smaller area. 

Almost all of our cores show strong 160-/im structure 
that is correlated with the 850-pim structure shown in the 
SCUBA data. Even those cores where SCUBA missed the 
emission peak show a surprising level of correlation, e.g. 
L1521C. The cores L1521A and L1709C were not detected 
with SCUBA, but the bounds of the observations are still 
shown in Figures [3] and |4] The 160-pm emission across 
those regions is relatively flat so it is possible that our non- 
detections were the result of there not being enough residual 
intensity to detect after the large scale emission had been 
chopped away. 

An exception to the 850/160/im pattern is the core 
L1696A (also known as Oph D). The small scale dumbbell- 
shaped SCUBA submillime tre emission closely m atches that 
seen in absorption at 7/im (Bacm ann et al.ll200Gt ) , but it ap- 
pears to be sat on a plateau of 160-/im emission that has 
no discernible substructure. Detailed radiative transfer cal- 
culations suggest that the southern part of th e core may be 
gravitationally bound ()Steinacker et allbOOSl ). 

L1521F, B68 and B133 all appear more isolated than 
the other cores, with strong central condensations whose 
closed contours extend almost to the edge of the maps. Half 
of the regions appear to be part of filamentary structures. 
Of the remaining images L1696A, L1498, and L1512 are on 
the edge of extended plateaux of emission while L1517 is a 
small cluster of star forming cores. L1517A and B are clearly 
seen at 160/im, but the SCUBA non-detection of L1517D is 
now seen not to be a distinct core at IBO/^m, so we discount 
it from our observations. 

Seven of our sources are associated with the L1521 
molecular filament that runs north-to-northwest from the 
west end of TMC-1 towards the L1495 region. A number of 
our fields touch each other so we have plotted their dis- 
tribution in Figure [5] The brightest of the L1521 cores, 
L1521F (also called MC 27), is separated from the main 
filament by half a degree and has been shown to be in an 
advanced evolu tionary state very close to the time of proto- 
star formation (Onishi. Mizuno. fc Fukuilll999l : rCrapsi et al.l 
[ioOSa V An embedded protostar has recently been reported 
in L1521F su ggesting that it has reached the end of its pre- 
stellar phase (|Bourke et al.|[2006l ) and is a very young, low- 
luminosity Clas s object somewhat similar to I RAM04191 
also in Taurus (|Andre. Motte. fc Bacmannlll999l ). 

L1521E is a triple-lobe d molecula r core situated at the 
end of the L1521 filament. iPaper v3 mapped the western 
lobe with SCUBA. This is the same lobe that has been 
identified as having a low chemical a ge based on com- 
parative studies of mole cular depletion l|Hirota et al.|[200^ : 
iTafalla fc Santiagdl2004l ). 

The two SCUBA peaks seen in the L1521SMM field 
seem to be part of a chain of sources that run east-west 
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Figure 3. Images of MIPS 160-/im (greyscale with black contours) and SCUBA 850-/jm (white contours) continuum emission towards 
the pre-stellar cores L1498, L1521 (C, SMM, B, A, D, F, and E), L1517, L1512, L1696A and L1709A. The smoothed black 160-Atm 
contours start at the level of the on-cloud background (see Table [2J and proceed in plus (solid lines) and minus (dashed lines) 2-cr 
increments. The lowest contour, the single dotted contour, marks the off-cloud background (see Tabic [2]l. The white 850-/im contours 
are at 45, 75, and 95 per cent of peak intensity. The mapped SCUBA area is bounded by a dashed white line. Only this bound is shown 
for SCUBA non-detections. 
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Figure 4. The pre-stellar cores L1689A, L1709C, L1689B, B68, and B133. Details as in Figure |3] 



across the map and are situ ated at the j unctio n between the 
clouds 16a and 16b in the iOnishi et al.fl|2002l ) H"C0+ sur- 
vey. The western-most SCUBA source is coincident with the 
position of a 2-/im source id entified as a brown dwarf, desig- 
nated J04210795+2702204 (|Guieu et al.ll2006t ). The eastern 
SCUBA source, which we deno te as L1521SMM, is one of 
our pre-stellar cores from iPaper VI The flux density from 
this source was measured in a 75-arcsec aperture to avoid 
confusion with the nearby object. 

The map of L1689A shows a filament running diagonally 
across the map with a break in it caused by a gap in the 
survey ed area. This is labelled Filament 1 bvlNutter et al.l 
(|2006t 1. There are two pre-stellar cores from lPaper Vil in this 
filame nt, but the L1689 SMM core lies under the coverage 
break. [Paper Vj onlv partially mapped L1689A due to time 
constraints, but other SCU BA observations have shown that 
the c ore is double peaked ijShirlev et al. I I2OOOI : iNutter et al.l 
I2OO6I ). We have added archival SCUBA data to our original 
map to produce the contours shown in Figure [3] 

L1689B is an elongated core which has been shown to 
be sharp-edged in the north-south direc tion based on mm 
emission and mid-IR absorption maps (| Andre et all 1 19961 : 
iBacmann et al.ll2000l ). The morphology seen in the Spitzer 
160-/^m image supports the conclusion drawn from mm and 
mid-IR data. 

B68 is a low-density circular core situated close to 
the Galactic Plane which has been identifi ed as a crit- 
ical Bonnor-Ebert sphere ()Alves et all l200lh . The bright 
source to the sou th-east is the starless core L55 SMM 1 
jVisser et al.ll200ll ). 



4 SPECTRAL ENERGY DISTRIBUTIONS 

We can add 160-/im fiux dens ity measu rements to our spec- 
tral energy distributions from lPaper Vl . However we are lim- 
ited by the extent of our SCUBA maps and cannot use the 
centroid flux densities presented in Tabled as they are not 
always coincident with the SCUBA positions. It is therefore 
necessary to measure flux densities in a common aperture 
centred on a common position for all wavelengths. 

Table [3] lists the 14 sources common to iPaper Vl . 

iPaper Vll and this paper, the positions of the SCUBA cores 
and their multiband flux densities measured in a 150 arcsec 
aperture centred on that position. Columns 6, 8 an d 9 list the 
ISOPHOT 70, 170, and 200- ^im flux densities from lPaper vl . 
table 4. Column 7 lists the Spitzer 160-/im flux density mea- 
sured at the SCUBA position (note this is a different po- 
sition and flux density than listed in Table (2}. Columns 10 
and 11 lis t the SCUBA 450 and 850-^m flux densities from 
IPaper v3 . tables 1 and 2. Four sources, L1521A, B & C, and 
L1709C were either undetected or missed by SCUBA. 

In order to enhance our short wavelength coverage, 
background subtracted fiux density measurements were 
made at the SED positions using 24- and 70- /xm Spitzer d&ta.. 
These were downloaded from the Spitzer archive. C2D legacy 
data were used for the Ophiuchus regions (R.A.=16''). Of 
the 14 pre-stellar cores with Spitzer data only 5 were de- 
tected at 70/xm and only L1521F was detected at 24^m. 
Columns 4 and 5 of Table |3] list the 24- and 70-/im flux den- 
sities and 1-a uncertainties or, where appropriate, 3-cr upper 
limits. 

At 24-^m several of the Ophiuchus cores were silhou- 
etted against brighter background emission. The 3-a upper- 



8 J.M. Kirk, D. Ward- Thompson and P. Andre 




30 28 26 24 22 4:20 18 



Right ascension (2000) 

Figure 5. The L1521 filament in Taurus. IRAS 100-/xm grayscale superposed with white MIPS 160-/im contours. The contour heights 
start at 40 MJy/sr and proceed at intervals of 15 MJy/sr thereafter. The dashed white line shows the extent of the MIPS data. The 
seven boxes show the subregions shown in Figure [3] that are included within this area. 



Table 3. SCUBA positions and multi-band flux densities with l-cr errors used in the calculation of the dust SEDs. The l-tr errors include 
contributions from statistical noise and instrumental calibration errors (20 per cent for MIPS, 30 per cent for ISOPHOT, 25 per cent 
for SCUB A 450/im, and 10 per cent for SCUBA 850/.tm). It was found necessary to revise a c ouple of ou r SCUBA flux densities from 
iPaper v3 . For example, slightly deeper 850-/im polarimetry observations of L1498 showed that iPaper v3 ha d under-estimate d the flux 
density at 850/im due to the orientation of the chop throw and the correction for this is included below (see iKirk et al.ll200(i) . 



Source Right Ascension Declination SED Flux Densities 

(2000) (2000) 24-/im ^Q-^lU^ gO-^tm 160-/im 170-Afm 200-^tm 450-Atm 850-Afm 

[mJy] [mjy] [Jy] [Jy] [Jy] [Jy] [Jy] [Jy] 



L1498 


04'i 


10"' 


52 


8^ 


-|-25°10'09" 


<1.3 


<30 


<1.0 


7.0±1.4 


10±3 


15±5 


13±3 


3.0±0.3 


L1521SMM 


04'' 


21'" 


00 


4" 


-|-27°02'32" 








0.86±0.20 








0.79±0.15 


L1521D 


04^ 


27m 


47 


0' 


-1-26° 17' 38" 


<3.3 


<61 




3.4±0.7 






14±4 


3.1±0.3 


L1521F 


04^ 


28™ 


39 


2^ 


-(-26°51'36" 


35±4 


290±60 




4.4±0.9 






12±4 


3.2±0.4 


L1521E 


04ft 


29"' 


13 


6= 


-|-26°14'05" 


<2.6 


<67 




0.66±0.23 






4.0±1.4 


1.4±0.2 


L1517A 


04^ 


55™ 


08 


9^ 


-l-30°33'42" 


<1.6 


46±18 


<1.2 


5.8±1.2 


9.6±2.9 


14±5 


4.8±1.3 


0.73±0.14 


L1517B 


04^ 


55™ 


17 


9^ 


-l-30°37'47" 


<1.6 


<46 


<0.78 


3.4±0.7 


5.9±1.8 


8.4±3.2 


12±3 


2.6±0.3 


L1512 


05'' 


04™ 


08 


4s 


-(-32°43'26" 


<1.4 


41±12 


<1.4 


9.9±2.0 


12±4 


16±5 


8.0±2.3 


1.4±0.2 


L1696A 


16'' 


28™ 


28 


9= 


-24°19'09" 


<13 


<590 


<2.1 


38±8 


39±12 


48±14 




6.3±0.6 


L1709A 


16'' 


30™ 


51 


9= 


-23°41'52" 


<4.4 


<250 


<1.6 


8.9±1.8 


10±4 


18±5 




1.5±0.2 


L1689A 


16'' 


32™ 


13 


2' 


-25°03'45" 


<4.0 


13000±3000 


35±10 


130±30 


130±40 


120±40 


19±6 


2.4±0.3 


L1689B 


16'' 


34™ 


48 


2' 


-24°38'04" 


<15 


<330 




35±7 


29±9 


29±9 


18±6 


3.1±0.4 


B68 


17'' 


22™ 


39 


2= 


-23°50'01" 


<2.3 


<75 


<2.8 


12±2 


12±4 


16±5 


<8.5 


1.3±0.2 


Bias 


19'' 


06™ 


08 


4= 


-06°52'22" 


<2.0 


110±30 


<4.5 


19±4 


23±7 


35±10 


16±5 


1.7±0.2 
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Table 4. SED fitting results. Columns 2, 3 and 4 list the tem- 
perature, peak wavelength, and luminosity (between 30/xm and 
3mm) of the best fit SED. Column 5 lists the reduced param- 
eter for the SED fit quoted to 1 s.f. It was assumed that (3=2 . 
Column 6 lists the equivalent SED temperatures from iPaper Vi 
L1521SMM has too few data points for a sensible determination 
of the parameter, but its data are consistent with the SED of 
a 9K greybody. L1521F has a second hotter dust component as 
shown in Figure (6] columns 3-5 list parameters for the colder of 
the two components while the parameter listed in Column 5 
is for the combined fit. 



Source 


T 


^peak 


LiR-MM 




T/so 




[K] 


[^tm] 


[L©] 




[K] 


L1498 


10.1±0.3 


290 


0.15 


0.1 


10 


L1521SMM 


g.oii.o 


320 


0.03 






L1521D 


9.0±0.3 


320 


0.11 


0.2 




L1521F 


9.3±0.3 


310 


0.12 


0.1 




L1521E 


8.1±0.4 


360 


0.031 


0.3 




L1517A 


12.6±0.4 


240 


0.093 


0.6 




L1517B 


9.3±0.3 


310 


0.10 


0.2 


10 


L1512 


12.1±0.3 


240 


0.14 


0.04 




L1696A 


11.6±0.4 


250 


0.46 


0.3 


10 


L1709A 


11.6±0.5 


250 


0.12 


0.2 




L1689A 


18.5±0.5 


160 


1.4 


0.3 


19 


L1689B 


12.5±0.5 


240 


0.32 


0.6 


11 


B68 


12.5±0.5 


240 


0.13 


0.04 


10 


B133 


12.8±0.3 


230 


0.49 


0.6 


13 



umn 3 of Table |4] lists the wavelength of peak emission, 
^peak, for each best fit SED. Column 4 lists the infra-red to 
mm luminosity, Lir-mm derived by integrating the best-fit 
SEDs between 30/im and 3mm. 

We have both ISOPHOT 170Atm and MIPS 160^m fiux 
densities for ten cores. All except one of the 160/im flux den- 
sities are either comparable with or slightly less than the 
170/im flux densities, in broad agreement with what would 
be expected for a 10-13K SED. The agreement in fiux den- 
sities for the Ophiuchus cores is particularly interesting as 
their peak intensities from Table [2] would indicate they are 
heavily into the MIPS saturated regime. 

The SED temperatures for L1521D and L15 21E are ap- 
proxim ately 9K, validating our assumption in iKirk et aP 
of thermal equilibrium between the dust temperature 
and the therm al rotational line temperatures measured by 
ICodella et al The temperatures for the three L1521 

cores are slightly lower than for the other cores. L1521F 
shows a short wavelength excess above the single cold SED 
fitted to the rest of the cores. A second hotter greybody can 
account for this excess as shown in figure |6] This second 
component has a temperature of 63±5K and a luminosity of 
~ 0.01 Lq making it one tenth as luminous as the extended 
core. This is consistent with the properties of the pro tostel- 
lar object detected in L1521F bv lBourke et al.l l|2006l '). 

The mean temperature of these prestellar cores is 12K 
and is significantly lower than the mean temperature (23K) 
of equivalent isolated cores with an embedded protostar 
iLaunhardt et al.lll997l '). 



limits derived from the 24-^m images varied according to the 
region and the scan speed. A typical 24-/im 3-a upper-limit 
for an isolated Taurus core observed using a medium-scan 
speed was -^l.GmJy (e.g. L1517B), whereas the 24-/^m 3-a 
upper-limit for some Ophiuchus cores observed using a fast- 
scan speed was an order of magnitude higher (e.g. L1689B). 

The data from Tableland selected data from lPaper Ilj 
are plotted on log-log plots as a function of wavelength in 
Figure [6] Fitted against the data from Table |3] and plotted 
in Figure |6] are a series of modified blackbody, or greybody 
curves. The monochromatic fiux density Si, of a greybody, 
at frequency u, radiated into solid angle Q is given by 



(1) 



where i3y,T is the Plank function, i/c is the frequency at 
which the optical depth is unity, Q is the solid angle of 
the aperture, / is the filling factor of the source within 
the aperture and /3 is the dust emissivity index. A grey- 
body was fitted through each set of data points by using x^- 
minimisation. I t was assumed that /3 = 2 and that v/uc = 1 
at 50 /xm, after IWard- Thompson et al.1 ()200 2fl. The parame- 
ters / and T were then mapped on a 1000 x 1000 grid of x^- 
Typical l-cr deviations from the best fit were found to be 
between 0.3 and 0.5K. 

Our best-fit parameters are listed in Tab le [4] along with 
the original ISOPHOT flt temperature from iPaper VI . The 
agreement between the two temperatures is seen to be good. 
The slight difference between the two is due to the addition 
of the 160-/im point which gives greater weight to the fit 
around the peak. There are no detected systematic differ- 
ences between the original and updated temperatures. Col- 



5 DISCUSSION 

5.1 Core Classification 

Recently there have been reports of very low lumi- 
nosity (L < 0.1 Lq) protostellar objects embedded 
within molecular co res that we r e prev io usly believed 
to b e pre-stellar jCrapsi et al.l l2005bl : iBourke et al.l 
I2OO6I ). These recently discovered protostars are simi- 
lar to some previous ly known Class protostars such 

as VLA 1623 jAndre. War d- Thompson, fc Barson- 

I993II, HH24MMS ll Ward- T hompson ct J] 



jsonyl 
1993 : 



Bontemps. Ward-Thornpson fc Andrei I199Q) . and IRAM 



04191 (|Dunham et al.1 120061 ). The core L1521F is now 
believed to harbour a young proto stellar object, based on 
its detection by Spitzer at 24/im l|Bourke et al.l [2OO6I ). A 
detection at this wavelength seems to be key in identifying 
embedded protostars in dense cores that were thought to 
be pre-stellar, due to the great sensitivity of Spitzer in the 
mid-infrared. 

We have confirmed the L1521F detection above and 
shown that this requires there to be hotter dust within 
L1521F - a signature of an internal heating source. We 
can determine from our combined Spitzer and SCUBA data 
whether any other pre-stellar cores contain warm dust, and 
therefore may also contain young protostars. 

Of th e 29 pre-s tellar cores detected by SCUBA at 
850-pim in iPaper Vll 23 had coincident MIPS obervations 
present in the Spitzer data archive, and of those 14 had us- 
able 160- /im data (as represented by the SEDs in Section |4]). 
One of the remaining cores, L1524, is actually protostellar 
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Figure 6. Spectral Energy Distributions for the cores listed in Table |3] Each plot, where available, shows the flux densities and l-cr 
calibration uncertainties or 3-(t upper-limits at 1.3mm (IRAM), 850 and 450-/im (SCUBA), 200, 170, and 90-Atm (ISOPHOT), and 160, 
70 and 24- /xm (MIPS). The solid curves are greybody fits to the data. A dashed curve indicates the presence of a second hotter dust 
component in the greybody fit. This is only seen in one core, L1521F, now believed to be protostellar, where the warmer component has 
T = 63 ± 5K. 



and had been misclassified. The remaining eight cores are 
L1686, L43, L234E, L63, L1148, L1155C, H and D. 

iKauffmann et al.l l|2005h identified a protostellar candi- 
date in L1148. However, their source, L1148-IRS, is not co- 
incident with our SCUBA detection and is a separate neigh- 
bouring object. We suggest the name L1148-SMM for the 



prestellar core to differentiate if from the protostellar core 
L1148-IRS. The separation between L1148-IRS and L1148- 
MMS is 200 arcsec, which equat es to ^0-3 pc g iven a dis- 
tance to L1148 of 325 pc (Str aizvs et all llQQal. Archival 
and C2D legacy 24-/im maps of the remaining 8 cores were 
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checked for the presence of embedded protostars, but no ob- 
vio us candidates were f ound . 

IStamatellos et aD ('2005') defined three criteria for dif- 
ferentiating young protostcllar cores from pre-stellar cores 
based on their SEDs (T > 15K, Xpeak < 170 /im) and 
submillimetre morphology (approximately circular within 
4000AU). One of our cores, L1689A, satisfies both condi- 
tions. This is also the core whose SED peaks furthest from 
the radio and which is most luminous. L1689A has two 
peaks in the submillimetre and is more c ondensed than the 
other cores. Based on the criteria listed bv lStamatellos et al.l 
l|2005l ). L1689A could have been considered a potential pro- 
tostellar candidate, but it was not detected by Spitzer at 
24^m, and so does not appear to contain any warm dust. 
Likewise, no other cores in our sample other than L1521F 
were detected at 24/^mwithin the cores' full-width at half- 
maximum (FWHM, see Table [Sjl at 160/im. 

Therefore, out of the 22 cores common to our Spitzer 
and SCUBA samples that were previously believed to be 
pre-stellar, only one contains a protostar. This appears to 
indicate that very low luminosity protostellar objects em- 
bedded within pre-ste llar cores are rare contrary to some 
early suggestions fe.g. lYoung et al.ll2004l '). 

We used the method of s tatist ical lifetimes 
l|Beichman et all Il986l : iLee fc Mve"r3 \W9^ to estimate 
the Hfetime of a prestellar core. In short, this method 
uses the relative abundance of sources to determine their 
relative statistical lifetimes in a general source population 
(the rarer an object the shorter its lifetime). One caveat 
to this method is that it requires a known lifetime for the 
entire system or of one part of it in order to calibrate the 
rest of the inferred lifetimes. For starless cores the reference 
lifetime of ~ 6 x 10^ years is inferre d from the lifetim e of 
the Class I protostellar phase (see iKirk et ahl |2005| . and 
references therein). 

From the starless core lifetime we estimated that the 
statistical lifetime for a prestella r core detectible in the sub- 
millimetre was ~ 3.0 X 10'' years (|Kirk et al.ll2005l ). This was 
based on a sample of 27 cores detected out of 50. A misclas- 
sification of only 1 out of 22 cores, as shown above, would 
shorten this estimate by ~5 per cent, which is well within 
the errors associated with th is lifetime (for a discussion see 
I Ward- Thompson et al.ll2006l ). Hence we see no reason to ad- 
just the pre-stellar core lifetime estimate based on the new 
Spitzer data. 

5.2 Mean Core Properties 

Table [5] shows the me an core pr operties derived from our 
SCUBA observations (jPaper Vj. table 4) and our Spitzer 
observations (this paper). The first line lists the sample size 
with the numbers in parenthesis referencing the 7 cores that 
are common to both samples (L1521F has been excluded 
from both). Lines 2-4 list the measured peak intensities and 
flux densities in a 150 arcsec aperture. A "bg" subscript 
denotes a background subtracted measurement or a quantity 
derived from that measurement. SCUBA observations are 
automatically counted as background subtracted because it 
uses a chopping secondary mirror to reject sky emission. 

Line 5 lists the mean of the core masses derived from 
the 150 arcsec aperture fluxes at the two wavelengths. When 
constructing an SED we assume that the same mass of gas is 



Table 5. Table of mean core statistics derived from tlie MIPS 
and SCUBA data. Column 2 lists the mean value derived from 
the Spitzer 160-/im observations for the thirteen prestellar cores 
discovered in this paper. Column 2 lists the equivalent mean pa- 
ram eter derived from the SCUBA 850-fim observations in Table 
4 of iKirk et all ||2005|) (except L1521F and L1582 whose greater 
distance skews the mean). The values in parenthesis are for the 7 
cores common to each sample. A "bg" subscript denotes a quan- 
tity that has had the background intensity subtracted. 



Spitzer WOfim SCUBA 850/im 



Sample size 


13 (7) 


10 (7) 






230 (170) 




MJy/sr 




86 (57) 


34 (27) 


MJy/sr 




22 (13) 


2.7 (2.2) 


Jy 


-^^ISOarcscc , bg 


1.4 (1.5) 


1.6 (1.3) 


Mq 


A^(H2)poak 


7.8 (5.4) 




10^2 cm 


-'V(H2)pcak,bg 


1.7 (1.4) 


4.5 (3.1) 


10^2 cm 


FWHMa 


0.125 (0.140) 


0.051 (0.058) 


pc 


FWHMb 


0.084 (0.086) 


0.028 (0.033) 


pc 


HWHMmean 


0.051 (0.054) 


0.016 (0.022) 


pc 


Aspect Ratio 


0.67 (0.69) 


0.59 (0.61) 






1.7 (1.7) 


0.5 (0.4) 


M© 


rBE 




~0.150 


pc 


^CC 


~0.120 




pc 



radiating at all wavelengths. If that assumption holds then 
we should expect to derive identical mass estimates from 
different points on the SED. That the two masses on this 
line are virtually identical is apparent proof of the afore- 
mentioned assumption. 

Lines 6 and 7 of Table[5]list the peak H2 column density 
before and after background subtraction. If the core is cen- 
trally condensed then we could expect the smaller 850-/im 
beam to preferentially sample a higher mean column density 
than the more extended 160-/im beam, which it does. 

Lines 8 and 9 list the mean deconvolved major and mi- 
nor full-width at half-maximum (FWHM) of the cores. For 
160/xm this is the half-maximum between the peak and the 
on-cloud background. Line 10 lists HWHMmean, the geo- 
metric mean of half of the FWHM ellipse and line 11 lists 
its aspect ratio. Line 12 lists MpwHM.bg the mean estimated 
mass enclosed by the FWHM ellipse. The Misoarcsec.tg and 
A^FWHM,6g are nearly identical at 160/xm because the mean 
core diameter translates to an angular size of ~150 arcsec 
at t he typical distance of these co res (140 parsec). 

IWard -Thompson et all (|2006l ) defined prestellar cores 
as "that subset of starless cores which are gravitationally 
bound and hence are expected to participate in the star 
formation process". F or the pres tellar cores observed with 
SCUBA M^ir ~ Mobs (iPaper Vj) and can therefore be con- 
sidered bound. Between the two wavelengths shown in Table 
[5]the observed radius and mass increase by a factor of three. 
For a virialised cloud M oc Av^R where Av is the turbu- 
lent linewidth of the gas. Therefor e, for a c onstant linewidth 
M is linear with respect to R (c.f. lLarsonlll98D . This is the 
same as the linear scaling seen between the cores as observed 
with SCUBA and as observed with Spitzer so it is reasonable 
to assume that the extra material seen by Spitzer is also in 
bound. 
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When comparing the FWHM size and mass of the cores 
we note that the mean SCUBA sizes and masses are in agree- 
ment with othe r SCUBA size/ma ss studies of nearby star 
formation (e.g. iMotte et al.l Il998l ) whilst the mean of the 
Spitzer sizes and masses are close to su bmillimetre studies 
of condensations in the Orio n region (e.g. iMotte et al.]l200ll : 
iNutter fc Ward-Thompsonl l2006) . 

With these observations and our SCUBA observations 
we are able to probe a spatial regime down to 2000 AU (the 
FWHM of the SCUBA 850^m at 140 parsecs). For the 13 
confirmed prestellar cores we predominantly see a one to 
one mapping of closed SCUBA contours onto closed Spitzer 
contours suggesting that we are detecting the same single 
coherent unfragmented object in both sets of observations. 
Most stars form in multiple systems (see iGoodwin et al.l 
l2006l . and references there in) with the implication that at 
some scale a single core must fragment into multiple proto- 
stellar fragments. That we see no such fragmentation over 
the scales that we ar e sensitive to supports the finding of 
ICoodwin et al.l (|2004l ) that fragmentation must occur below 
the 2000 AU scale and/or at a later core evolutionary stage 
(such as the Class stage). 

Fitting of Bonner-Ebert (BE) spheres to SCUBA maps 
gives an implied pressure bounded outer r adius oirsE ■= 3 x 
lO'^ AU (0 150 p arcecs; line 13 of Tabled (|Evans et al.ll200ll : 
iKirk et ahlbOOSh . This is approximately equal to the major 
FWHM of the Spitzer cores and the radius of dense cores 
mapped in C^*0 jMvers et al.lll99ll : IOnishi et al.lll998l '). The 
mean HWHM of cores in the Spitzer observations is ap- 
proximately equal to the size of prestellar c ores measured 
from N2H-f (0 054 parsecs: ICaselh et al.ll2002f) a r id H^ ^CO-f 
(0.048 parsec; lOnishi et al.ll2002t ). lYoung et all (|2006f l used 



the SCUBA scan mapping mode to map larger areas over 
a selection of prestellar and protostellar cores. Their maps 
showed that the area of emission detected from cores with 
SCUBA was generally smaller than 2 arcmin across (the 
width of the SCUBA field of view). 

From the above radii we can define a "characteristic" 
radius for the full extent of a prestellar core. This is ~0.15 
parsec and is the BE pressure bounded radius and the radius 
as measured from CO observations. At the distance to these 
sources this radius is equivalent to ~ 3.5 arcmin. Given the 
complexity of the lowest on-cloud contour in Figures [3] and 
|4] it is nearly impossible to measure a full extent radius. 
However, by inspection, the regions typically have a closed 
contour radius, Tcc, of ~3 arcmin (0.12 parsecs; line 14 of 
Table [5]), comparable to the 3.5 arcmin BE radius. Then, 
within this bounding radius there is a centrally condensed 
region of radius ~ 0.05 parsec as delineated by observations 
of high density gas tracers and the FWHM of far-IR dust 
emission. 



6 CONCLUSION 

In this paper we have presented Spitzer archive data of a 
sample of 18 molecular cloud cores that were believed to 
be pre-stellar in nature. The sample was selected from the 
52 cores that we had previously observed with SCUBA. We 
have produced new images of the cores at 160/im and we 
have measured their flux densities at the other MIPS wave- 
lengths of 24 and 70fim to produce full SEDs for the cores. 



when combined with our previous data. Comparison of back- 
ground intensities and measured flux densities showed con- 
sistency between the absolute and relative MIPS ITO/^mand 
ISOPHOT 160/imcalibration. Four of the cores were non- 
detections in our SCUBA sample. We refined our previous 
fits to the remaining 14 core SEDs and produced tempera- 
ture estimates for the cores. Most lie in the range 10-13K. 
Only one core was detected at 24/im and seen to need an 
additional temperature component to fit its SED - L1521F. 
This core has recently been reported to contain a protostar, 
and we confirm this assertion. However, no other core was 
detected at 24/im and no other core required a hotter dust 
component to fit its SED. Hence we deduce that the number 
of misclassified pre-stellar cores is probably small, and the 
resultant change in the estimated lifetime of the pre-stellar 
core stage is <5 per cent. 
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